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' rections. The results clearly favor distribution amplitudes close to the asymp- 
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tion on the elastic form factor of the pion and on its valence quark distribution 
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photon decay widths, we determine the 77 and the 7/ decay constants and the 
77-7/ mixing angle. 
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I. INTRODUCTION 



Hadronic form factors at large momentum transfer provide information on the con- 
stituents the hadrons are made of and on the dynamics controlling their interactions. There- 
fore, the form factors always found much interest and many papers, both theoretical and 
experimental ones, are devoted to them. Recently a new perturbative approach has been 
proposed by Botts, Li and Sterman which allows to calculate the large momentum 

behavior of form factors. In this new approach, which one may term the modified hard 
scattering approach (HSA), the transverse degrees of freedom as well as gluonic radiative 
corrections — condensed in a Sudakov factor — are taken into account in contrast to the stan- 
dard perturbative approach [|J. Applications of the modified HSA to the pion's and nucleon's 
form factor p],|3|,|5||| revealed that the perturbative contributions to these form factors can 
reliably and self-consistently (in the sense that the bulk of the contributions is accumu- 
lated in regions where the strong coupling a s is sufficiently small) be calculated. It turned 
out, however, that the perturbative contributions are too small as compared with the data. 
Responsible for that discrepancy might be omitted contributions from higher order pertur- 
bation theory and/or from higher Fock states. Another and perhaps the most important 
source of additional contributions to the form factors are the overlap of the soft wave func- 
tions for the initial and final state hadrons ]7|]. All these contributions are inherent to the 
standard as well as to the modified HSA but, with very few exceptions, have not been con- 
sidered as yet. 

In this paper we are going to apply the modified HSA to pseudoscalar meson-photon tran- 
sition form factors for which data in the few GeV region is now available fSHlQfl. These 
transition form factors are exceptional cases since there is no overlap contribution in con- 
trast to the elastic hadronic form factors and, to lowest order, they are QED processes. 
QCD should only provide corrections of the order of 10 — 20% and higher Fock state con- 
tributions are suppressed by additional powers of a s /Q 2 . For these reasons one may expect 
the modified HSA to work for momentum transfer Q larger than about 1 GeV. The low Q 
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behavior of the transition form factors has been investigated in a Salpeter model in [|TT|1 . 
Input to calculations within the modified HSA are the hadronic wave functions which con- 
tain the long-distance physics and are not calculable at present. However, the pion wave 
function required for the calculation of the 71-7 transition form factor, F 7r7 (<5 2 ), is rather 
well constrained. A reliable quantitative estimate of the ir-'j transition form factor can 
therefore be made, as we will discuss in Sec. |J. We extend this calculation in Sec. [I I J to 
the 77-7 and 7/-7 form factors and determine the decay constants and the mixing angle for 
pseudoscalar mesons. Sec. fV| is devoted to a discussion of the elastic form factor of the 
pion and to its valence quark distribution function. It will turn out that the same pion 
wave function as is used in the calculation of the 7r-7 transition form factor, also leads to a 
reasonable description of the valence quark distribution function and the elastic pion form 
factor provided the overlap contribution is also taken into account in the latter case. We 
recall that the Drell-Yan-West overlap formula is the starting point of the derivation of 
the hard scattering formula for the elastic pion form factor. The paper terminates with a 
few concluding remarks (Sec. |Vp. 



II. THE 7T-7 TRANSITION FORM FACTOR 

Adapting the modified HSA to the case of ir-'y transitions we write the corresponding 
form factor as 

r d 2 b -> - 

Fn~,(Q 2 ) = J dxi-^— ^0(^1,-^1 ) T H {xi,b u Q) exp[-S(x u bi,Q)] , (2.1) 

where 61 is the quark-antiquark separation in the transverse configuration space. x\ and 
x 2 = 1 — X\ denote the usual momentum fractions the quark and the antiquark carry, 
respectively. Tjj is the Fourier transform of the momentum space hard scattering amplitude 
to be calculated from the Feynman diagrams shown in Fig. [|. Neglecting the quark masses 
and the mass of the pion, the hard scattering amplitude T# in momentum space reads 



where the charge factor C n is (e u 2 — e^ 2 ) / a/2, e, denotes the charge of quark i in units of 
the elementary charge. The Fourier transform of this amplitude reads 

T H (xi, h, Q) = ^^K (^Qbi), (2.3) 

7T 



where K is the modified Bessel function of order zero. Strictly speaking, (|2.3|) represents 



twice the Fourier transform of the first term in (2.2). Because of the symmetry of the pion 
wave function under the replacements X\ «-> X2, which appears as a consequence of charge 
conjugation invariance, the second term leads to the same contribution as the first term after 



integration over the transverse separation b%. The Sudakov exponent S in ( |2.1| ) comprising 
the gluonic radiative corrections, is given by 

Six,, Q, h) = s(x u Q, h) + s(x 2 , Q, b{) - A In ^^BL (2.4) 

where a Sudakov function s appears for each quark line entering the hard scattering am- 
plitude. The last term in (|2.4|) arises from the application of the renormalization group 
equation (/3 = 11 — 2/3n/). A value of 200 MeV for Aqcd is used throughout and t is taken 
to be the largest mass scale appearing in T H , i. e., t = max( v /x^Q, l/&i). For small b± there 
is no suppression from the Sudakov factor; as b\ increases the Sudakov factor decreases, 
reaching zero at 61 = 1/Aqcd- For even larger b\ the Sudakov factor is set to zero. The 
Sudakov function s is explicitly given in |1].|2[ . 

The quantity \l/o appearing in (|2.1| ) represents the soft part of the transverse configuration 
space pion wave function, i. e., the full wave function with the perturbative tail removed 
from it. Following |J we write the wave function as 

* (zi, h) = ^= (j){x x ) E(V^fci)- (2.5) 

It is subject to the auxiliary conditions 

E(0) = 4vr, f 1 dxt 0(xi) = 1. (2.6) 

J 

The wave function does not factorize in x\ and bi, but in accord with the basic properties 
of the HSA ||T^|,|IH &i-dependence rather appears in the combination *Jx\X2 b±. The 



transverse part of the wave function is assumed to be a simple Gaussian 



fl(y/xiX2 bx) = Air exp Xix 2 b\/Ac 



(2.7) 



In it is shown how (77?) is related to the equal time harmonic oscillator wave function. 
More complicated forms than (|2.7|) (e. g., a two- humped shape of the momentum space 
wave function) are proposed in [113 on the basis of dispersion relations and duality. At large 
transverse momentum, however, the soft momentum space wave function should behave like 



a Gaussian ||13|| . The examination of a number of examples corroborates our expectation 



that forms of E other than ( [2.71) will not change the results and the conclusions presented 
in our paper markedly. 

The two free parameters contained in our ansatz, namely the value of the wave function at 
the origin and the parameter a controlling the transverse size, are well fixed by the decay 
processes 7r + — > /U + z/ M and tt° — > 77 providing the relations |12| 



dxi ^oO^iA = 0) - n 
V6 



dx\ d b\ ^(^l) b\) 



(2.8) 



where / 7r (= 130.7 Me V) is the usual pion decay constant. The first relation is automatically 
satisfied by our ansatz ( |2.5[ ), whereas the second relation is used to fix the parameter a. For 
the distribution amplitude, 0, we use the asymptotic form 



<f>As{xi) = 6x1X2 



(2.9) 



and alternatively, as a representative of strongly end-point concentrated distribution ampli- 
tudes, a form proposed by Chernyak and Zhitnitsky fl4 



hzixi) = 30xia;2 (xi - x 2 ) 



(2.10) 



It can be shown (see 0] and references therein) that the distribution amplitudes are subject 
to evolution and can be expanded over Gegenbauer polynomials which are the eigenf unctions 
of the evolution equation for mesons 



0(xi) = <pAs{xi, 



n=2 



tt', 



a 



:(Mo) / 



X 2 



(2.11) 
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fiF is a scale of order Q at which soft and hard physics factorizes and /j,q is a typical hadronic 
scale of order 1 GeV. Charge conjugation invariance requires the odd n expansion coeffi- 
cients B n to vanish. Since the 7„ are positive fractional numbers any distribution amplitude 
evolves into 4>as{%) asymptotically, i. e., for ln(Q//z ) — > oo. The asymptotic distribution 
amplitude itself shows no evolution. In the representation (2.11) the CZ distribution ampli- 



tude is given by B2 = 2/3 and B n = for n > 2. Its evolution can safely be ignored in the 
very limited range of momentum transfer we are interested in. 

In the standard HSA the distribution amplitude ( [2.10|) leads to a prediction for the elastic 



pion form factor in apparent agreement with experiment at the expense, however, of the 
dominance of contributions from the end-point regions, X\ — > or 1, where the use of the 



perturbative QCD is unjustified as has been pointed out by several authors [pL5|,|T6|| . It is 
controversial whether or not (|2.10|) is supported by QCD sum rules. We do not want to 
enter that discussion but consider ( p.lOj ) as an example whose significance is given by its 
frequent use. 



From the ir° — > 77 constraint (|2.8p we find for the transverse size parameter a the values 



861 MeV for the asymptotic wave function and 673 MeV for the CZ wave function. These 
values for a correspond to the following characteristic properties of the pion's valence Fock 
state: 0.21 (0.32) fm for the radius, 367 (350) MeV for the root mean square transverse mo- 
mentum and 0.25(0.32) for the probability when the asymptotic (CZ) wave function is used. 
The small radius of the valence Fock state has a lot of implications in hard processes [|17j . 



Numerical results for the transition form factor F n ^(Q 2 ) obtained from fl2.1|) are displayed in 
Fig. 0. We emphasize that there is no free parameter to be adjusted once the wave function 
is chosen. Obviously the results obtained from the asymptotic wave function are in very 
good agreement with the CELLO data [|1(]] as is also indicated by a x 2 value of 6.2 for the 
five data points; there is not much room left for contributions from higher order perturba- 
tive QCD and/or from higher Fock states. The results obtained from the CZ wave function 
overshoot the data significantly (x 2 = 44.5). Thus the comparison of both the results with 
the data on the 7r-7 transition form factor evidently ends in favor of the asymptotic wave 



function. Of course, the limited quality and quantity of the data allows mild modifications of 
the asymptotic wave functions without worsening the agreement between theory and exper- 
iment considerably. For example, if we follow Brodsky et al. [I^j and multiply ( |2.9| ) by the 



exponential exp(— a 2 m 2 q / X1X2) where the parameter m q represents a constituent quark mass 
of, say, 330 MeV, we find similarly good results from this modified asymptotic wave function 
(termed BHL subsequently) as from (|2.9|) itself. On the other hand, strongly end-point 
concentrated wave functions are clearly disfavored. Admixtures of the second Gegenbauer 
polynomial enhancing the end-point regions, must be small; already a value of 0.1 for the 
expansion coefficient B2 deteriorates the agreement with the experimental data substantially 

(x 2 = 10). 

The standard HSA @,|18[ predicts for the 7T-7 transition form factor 

F^{Q 2 ) = V2/3 Mx^ 1 } Q- 2 . (2.12) 

The bracket term denotes the xf 1 moment of the distribution amplitude. This moment 
receives the value 3 and 5 for the asymptotic and the CZ distribution amplitude respectively. 
Obviously, the standard HSA, while exact at large Q, fails to describe the data in the 
few GeV region. It does not provide the substantial Q-dependence the data for Q 2 F TJ 
exhibits. This is to be contrasted with the modified HSA in which the QCD corrections, 
condensed in the Sudakov factor, and the transverse degrees of freedom provide the required 
Q-dependence. Asymptotically, i. e., for \n(Q/fi ) — > 00, the Sudakov factor damps any 
contribution except those from configurations with small quark- ant iquark separation and, 
as the limiting behavior, the QCD prediction [[18l , |T9|l F nry — > V2f n Q~ 2 emerges. 
The Q-dependence of the form factor can be parameterized as 

F 7T1 (Q 2 )=A/(l + Q 2 /s ) (2.13) 

where for the constant Sq a value of about 0.6 — 0.7 GeV 2 is required by the data. In vec- 
tor meson dominance models so is to be identified with the square of the p-meson mass 
(0.59 GeV 2 ). Brodsky and Lepage P§ propose the formula ( |2.13D as an interpolation be- 



tween the two limits, F 7rJ (Q 2 =0) = (2y/2 tc 2 /^) 1 known from current algebra and the above 



mentioned asymptotic behavior \[2f 1 ,Q~ 2 , hence s = ^ 2 f^ = 0.67 GeV 2 . We stress that 
this interpolation formula is not derived within the standard HSA. In QCD sum rule analy- 
ses so is related to the pion's duality interval POfl . O ur results respect ( |2.13|) approximately 



and so represents the net effect of the suppressions due to both the Sudakov factor and the 
&i-dependence of the wave function and of T#. 

We close the discussion of the 7r-7 transition form factor with a reference to a paper by 



Szczepaniak and Williams plf in which soft, nonperturbative corrections to the standard 
HSA are estimated. These corrections bear resemblance to the intrinsic transverse separation 
effects considered by us. 

III. THE 77-7 AND r/- 7 TRANSITION FORM FACTORS 



We are now going to generalize ( [2.1|) to the cases of 77-7 and rf-'y transitions. We start 



with the SU(3) basis states, r] 8 and 771, and employ the usual mixing scheme 

| 77) = cos iip I r] 8 ) — sin iip I r/i) 

I 77') = sin lip I 77 8 ) + cos-ip | 771). (3.1) 

Insertion of this scheme into the 77-7 and rj'-j matrix elements of the electromagnetic current 
leads to relations between the physical transition form factors and the 7/3-7 and the 771-7 
ones. The latter form factors can be calculated analogously to the tt-j case. For the r/ 8 and 
77! wave functions we use the same ansatz as for the pion (i = 1,8) 

fi 

61) = ^= 0i(xi) T,i(y/xix 2 bi) (3.2) 

and assume that, except of the decay constants fi, the two wave functions are identical to 
the asymptotic pion wave function, ( |2.7| ) and (|2.9|) . In the hard scattering amplitude ( |2.2| ) 
the charge factor of the pion is to be replaced by either 

C 8 = (e u 2 + e d 2 - 2e s 2 ) /V6 (3.3) 

or 



G x = (e u 2 + e d 2 + e 2 ) /Vs. 



(3.4) 



Furthermore we correct for the rather large masses of the ^-mesons (mg = 566 MeV, mi = 
947 MeV, see [0]). Hence, in the transverse separation space, the hard scattering amplitude 
reads 



f H {xi,bi,Q) = — K (J x 2 Q 2 + xxx 2 mi 2 bi). 
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(3.5) 



Inserting ( [3.2|) and (|3.5|) into (2T) we can compute the F mi form factors and, then, using the 
mixing scheme ( |3.1|) , the 77-7 and rf-j transition form factors. However, we also need for this 
calculation the values of the decay constants and the mixing angle. Since these parameters 
are not known with sufficient accuracy we will change our attitude and, encouraged by the 
success of the modified HSA in the 71--7 case, try to determine them. Admittedly, additional 
information is required for this task since the 77-7 and rj'-'y transition form factors do not 
suffice to fix the three parameters; for any value of the mixing angle a reasonable fit to the 
data is obtained. The necessary extra information is provided by the two-photon decays of 
the 77 and 7/. Adapting the PCAC result for the tt° — > 77 decay to the 77 and rj' case with 
proper mixing at the amplitude level, one finds for the decay widths 

9a 2 



r(»7 -» 77) 

r(r/ -> 77) 



16^™" 
9a 2 



16^ v 



—— cos Vp — — sin up 

J8 Jl 

Cg . Q . C\ 

—— smvp + — cosi/p 

J8 Jl 



1 2 



(3.6) 



The experimental values for the decay widths are p3|: 



r(rj -> 77) = 0.510 ± 0.026 keV 



r(V -> 77) = 4.26 ± 0.19 keV 



(3.7) 



For obvious reasons we have only quoted the PDG average of the two-photon measurements 
of the 77 — > 77 width. The value of 0.324±0.046 keV obtained from the Primakoff production 



measurement will not be used in our analysis. 

We evaluate the three parameters, fx, /g and $p, through a combined least square fit to 
the data on the form factors and the decay widths. The parameters acquire the following 
values: 
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fx = 145 ± 3 MeV, f 8 = 136 ± 10 MeV, $ P = -18° ± 2° (3.8) 

and the x 2 value is 14.8 for the 18 data points. The values of the parameters and their 
errors are only correct provided the r] and rj' wave functions are at least approximately given 
by ( |3.2| ), ( |2.7| ) and (O). Since f\ and fs have rather similar values nonet symmetry of the 



wave functions holds approximatively. The decay constants of the physical mesons are: 

f v = 175 ± 10 MeV; f v , = 95 ± 6 MeV. (3.9) 

The quality of the fit can be judged from Fig. ^| where fit and data |8|-[T0| for the transition 
form factors are shown. As expected from the very good value of the \ 2 the agreement 
between theory and experiment is excellent. The computed values for the decay widths are 
r (v -> 77) = 0.5 keV and T (r/ -» 77) = 4.17keV. 

Our value for the mixing angle is compatibel with other results. For instance, a value of 
—20° ± 4° is obtained from chiral perturbation theory ||22|| . An analysis of two-meson final 
states produced in proton-antiproton annihilations seems to provide a similar value for the 
mixing angle (— 17.3°±1.8°) [25|. Gilman and Kauffman |26|] found from a phenomenological 



analysis of various decay processes and of ir~p scattering that a value of about —20° is 
favored. From a similar analysis but under inclusion of constituent quark mass effects 
Bramon and Scadron |27j] obtained ftp = —14° ± 2°. 

From chiral perturbation theory Gasser and Leutwyler |22] predicted a value of 170 ±7 MeV 
for the t? 8 decay constant whereas Donoghue et al. found 163 MeV. While both the values 



are compatible within the uncertainties of the chiral perturbation theory, they are larger than 
our result. In order to see whether or not such large value is definitively excluded in our 
approach we repeat the combined fit, keeping / 8 at the value of 163 MeV. The resulting fit 
is not as good as the precedent fit but still of acceptable quality, the value of x 2 is 20.7. It 
provides: fx = 143±3MeV, $ P = -21°±1° and hence /„ = 203±2MeV, /„/ = 76±5MeV. 
The results for the form factors are almost as good as before. The resulting decay widths 
are r (77 — ^ 77) = 0.47 keV and T (7/ — > 77) = 4.20 keV, i. e., imposing larger values upon 
fs forces the two-photon decay width of the 77 towards smaller values closer to the result 
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of the Primakoff experiment ||24j| . In view of the experimental uncertainties in the 77 — > 77 
and of the moderate difference in the \ 2 we cannot exclude the possibility of a as large 
as 163 MeV although a value around 140 MeV is favored from our analysis. In contrast to 
fs the other two parameters, f\ and $p, are tightly constrained. We stress that the quoted 
errors are only those obtained in the statistical analysis. By no means they reflect the full 
uncertainties of the parameters which are rather represented by the differences between the 
two sets of parameters. In |9|Jl0t| the rj and rj' decay constants have been determined from 
fitting pole formula? analogous to (|2.13|) to the transition form factor data. The values 



obtained for the decay constants differ from our ones considerably. Mixing is not taken into 
account in these analyses. 

It is often speculated upon a gluon admixture to the rji 



m = y U u + dd + ss + eG\ /V3 + e 2 . (3. 10) 

Allowing for that component which only changes the singlet charge factor (|3.4j) through the 
normalization, and repeating the fits to the data, we do not find any evidence for a sizeable 
gluon admixture to the r/i (e = ± 0.2). 

Finally we would like to comment on another mixing scheme frequently discussed in the 



literature (see, for instance, p6] , |27| , p9f ) , in which one uses quark states, uu + dd and ss, as 
the basis. Along the same lines as discussed above one may calculate the transition form 
factors using this mixing scheme. However, since we do not know well enough the masses 
and the wave functions of the quark basis states we refrain from carrying out that analysis. 



IV. DO WE KNOW THE PION'S WAVE FUNCTION? 

In this section we are going to put together the available information upon the pion's 
wave function. For this purpose we call to mind previous analyses of the elastic form factor 
of the pion and of its valence quark distribution function and discuss the results of these 
analyses in the light of our observations made on the 7T-7 transition form factor. 
Within the modified HSA the pion form factor is to be calculated from the relation [@,|]|| 
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F/ ert (Q 2 ) = f Q ^0- J°°J% d% %( yi , b 2 ) f H {x x , yi, Q, K h, t) ^o(xi, -h) 

x exp [S(x!, Q, bi, t) - S( Vl , Q, b 2 , t)] (4.1) 

where the hard scattering amplitude Tg is given by the product of two zeroth order Bessel 
functions arising from the gluon and the quark propagators in the elementary Feynman 
diagrams 

f H {x uVl ,QMM) = Aas ^ Cp K (^xmQb2) Koi^Qlh + hl) (4.2) 

7T 

Cp{= 4/3) is the color factor and the Sudakov exponents S are defined in ( |2.4| ). t, 
the largest mass scale appearing in the hard scattering amplitude, is now given by 



max( v /xi?/iQ, 1/bi, l/b 2 ). Equation (|4.1|) implies two angle integrations, say the integra- 
tion over the direction of b 2 which simply provides a factor of 2n, and the integration over 
the relative angle between b\ and b 2 . The latter integration can be carried out analytically 
by means of Graf's theorem 

^jdtp K (c \h + b 2 \) = Oih - 6 2 )^o(c6i)/ (c6 2 ) + 6(b 2 - b^K^cb^ch), (4.3) 

where 6 is the usual step function and Iq{x) = Jo(ix); Jq is the Bessel function of order zero. 
Thus a four dimensional integration remains to be carried out numerically. 
The perturbative contributions to the elastic form factor obtained from both the wave func- 
tions, the asymptotic and the CZ one, are compared to the data |5[J in Fig. [|. Both the 
predictions, while theoretically self-consistent for momentum transfers larger than about 
2 GeV, fall short of the data. One may wonder about the reason for the failure of the mod- 
ified HSA in the case of the elastic form factor while it works so well for the n-j transition 
form factor. At this point we remind the reader that the elastic form factor also gets a 
contribution from the overlap of the initial and final state soft wave functions $o & A 
contribution of this type does not exist for the transition form factor. Formally the pertur- 
bative contribution to the elastic form factor represents the overlap of the large transverse 
momentum (k±) tails of the wave functions while the overlap of the soft parts of the wave 
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functions is assumed to be negligible small at large Q. Since the soft wave functions are 
explicitly used in applications of the modified HSA it is natural, even obligatory, to examine 
the validity of that presumption. If the overlap of the soft wave functions turns out to be of 
substantial magnitude it must be taken into account in the calculation of the elastic form 
factor for consistency! According to |7| the overlap of the soft transverse configuration space 
wave functions reads 

F^°f\Q 2 ) = ±-Jdx l j d\exp [-(1 - x 1 )b 1 ■ q ± ] |* (*i, &i)| 2 , (4-4) 



where Q 2 = qj_. As the inspection of (|4.4j ) reveals only the end-point region 1 > 
x\ > 1 — 1/ 1 (Q(b\) l l 2 ) contributes at large Q, i. e., only configurations where the photon 
interacts with a parton carrying almost the entire momentum of the pion. Higher Fock 
states provide similar overlap terms. Obviously with an increasing number of partons shar- 
ing the pion's momentum it becomes less likely that one parton carries the full momentum 
of the pion. Therefore higher Fock state overlap terms are strongly suppressed at large Q. 
Consider for example the qqg Fock state and assume that the corresponding wave function 
is a suitable 3-particle generalization of Q2.5|) , (2.7) and ( |2.9|) where the asymptotic form 



of the qqg distribution amplitude is ~ £1X2X3 (^3 refers to the gluon) |3T|. It is easy to 
convince oneself that the qqg overlap term behaves as 1/Q 12 at large Q; above Q — 2 GeV 
it can be neglected to any degree of accuracy. 

Soft contributions calculated from (|4.4j) are shown in Fig. [|. The results obtained from the 
asymptotic wave function obviously have the right magnitude to fill the gap between the 
corresponding perturbative results and the data.[] The broad flat maximum of the overlap 
contribution simulates a Q-dependence of the form factor which, in the few GeV region, 
resembles the (^"dependence predicted by the dimensional counting rules. For the CZ wave 
function the overlap contribution exceeds the data significantly; the maximum value of Q 2 F n 



1 We refrain from showing the sums of soft and perturbative contributions because some double- 
counting might be implied. 
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amounts to 2.1 GeV 2 and is located at Q = 7.4 GeV. We consider this result as a serious 
failure of the CZ wave function. 

At large Q the overlap contributions from our wave functions are suppressed by 1/Q 2 as 
compared to the perturbative contribution. At which value of momentum transfer the tran- 
sition from the dominance of soft to hard contributions actually happens depends on the 
end-point behavior of the wave function sensitively. Already a moderate modification of 
the wave function may shift that value considerably. For the asymptotic wave function the 
transition takes place at the rather large value of about 10 GeV while for the BHL wave 
function the additional exponential effectuates a sharper decrease of the soft contribution 
and the hard contribution takes the control at about 5 GeV (see Fig. f|). However, below 
3 GeV the contributions from both the wave functions do not differ much. They also provide 
similar results for the n-j form factor as we already mentioned. 



Large overlap contributions have also been observed by other authors [p| p!5| , p2] , p3| . Thus 
the small size of the perturbative contribution finds a comforting explanation a fact which 
has already been pointed out by Isgur and Llewellyn-Smith | 15| . Were that contribution 



(including higher order perturbative corrections) much larger the existing large overlap con- 
tributions would be inexplicable. This conclusion seems to be compatibel with calculations 
of the one- loop corrections to the perturbative contribution [0,0 • 

The structure function of the pion or rather its parton distribution functions q v {x) offer 



another possibility to test our wave functions against data. As has been discussed in |12 



and p6f the parton distribution functions are determined by the Fock state wave functions. 



Each Fock state contributes through the modulus squared of its wave function integrated 
over transverse momenta up to Q and over all fractions xi except those pertaining to par- 
tons of the type v. Obviously the valence Fock state wave function only feeds the valence 
quark distribution function q\ = q\. Since each Fock state contributes to the distribution 
functions positively, the inequality 



y i 





Q d 2 k ± 


j dx\ J 


167T 3 



V(xi,k ± )\ 2 S(xt-x) (4.5) 
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holds. In order to calculate the Q dependence of the valence quark distribution function 
we need to know the full pion wave function, its perturbative tail is responsible for the 
evolution behavior. Yet the soft wave function ( |2.5| ) provides the bulk of the valence Fock 
state contribution to the distribution function and, consequently, in the few GeV region this 
contribution should respect the inequality (f4.5p. Hence 

Since, on the average, higher Fock state partons possess smaller values of x than partons 
from the valence Fock state the inequality should become an equality for x — > 1. For other 
values of x many other Fock states contribute to the distribution function in general, a fact 
which has to be contrasted with exclusive reactions to which only the valence Fock state 
contributes at large Q. 

In Fig. [5] the results for the valence quark structure function obtained from the asymptotic 
and from the CZ wave functions are shown and compared to the parameterization 

xq^(x) = Ax Sl (l -xf 2 . (4.7) 

The constant A is determined by the requirement Jq q\ (x) dx = 1. Sutton et al. analysed 
the recent data on the Drell-Yan process ir~p — > fi + fi~X and determined the values of the 
other two parameters; they quote 8\ = 0.64±0.03 and (the averaged value) 5 2 = 1.11 ±0.06 
at Q = 2 GeV. In Fig. |5] the parameterization of |37j is displayed as a band indicating 
the uncertainties of it. The width of the band does not take into account the uncertainties 
induced by theoretical assumptions underlying that analysis (e. g., a. K factor). In particular 
near x = 1 the width of the band may be underestimated since the parameterization of the 
proton structure function, used as input in that analysis, is an extrapolation for x > 0.75, 
i. e., it is not supported by data. 

The comparison of the parameterization ( [4.7| ) with the predictions obtained from our wave 
functions clearly reveals that the asymptotic wave function respects the inequality ( [4.51 ). 
The little excess around x = 0.9 is perhaps a consequence of the above mentioned extra 
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uncertainties in the parameterization ( 4/7) and/or of small deviations of the wave function 
from the asymptotic form. On the other hand, the CZ wave function exceeds the result of 
f3?|j dramatically at large x. The above-mentioned uncertainties can not account for that. 
Thus, the CZ wave function fails again and is, therefore, to be rejected. This conclusion has 
already been drawn by Huang et al. [3B|. However, our analysis of the 71-7 transition form 
factor, a quantity which has not been considered in []3q| , strengthens the evidence for the 
asymptotic wave function and against the CZ one. 

Frequently the pion form factor is calculated via the overlap formula ft4.4|) using a wave 



function normalized to unity (e. g., |3^j3^|). It is customarily asserted that such a wave 
function describes the binding of constituent quarks in a pion. While such a wave function 
may provide an overlap contribution to the pion form factor in fair agreement with the 
data at all Q 2 and may also respect the 7r — > fiv constraint it is hard to see how a conflict 
with the data for the other three reactions we are considering, can be avoided. In order 
to illustrate the difficulties arising from a wave function normalized to unity, let us look to 
the asymptotic wave function again. Duplication of the parameter a normalizes the wave 
function to unity and then the following results are obtained: the contribution to the valence 
quark distribution function is four times bigger than that one shown in Fig. [|, the 7T-7 form 
factor becomes too large (x 2 = 49) and the 7r — > 77 constraint is badly violated. The pion 
form factor on the other hand is well described. We finally note that the BHL wave function 
(normalized to unity) leads to results for the elastic form factor which are very similar to 
those presented in []32||. Yet the deficiencies in the other reactions remain. 



V. CONCLUSIONS 

We summarize our findings about reactions involving pions as the only hadrons: The 
asymptotic wave function as the only phenomenological input leads to a successful descrip- 
tion of four exclusive processes and is compatible with the band on the pion's valence quark 
distribution function. The four processes are, on the one side, the leptonic decay of the 
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charged pion and the two-photon decay of the uncharged pion fixing the parameters of the 
wave function and, on the other side, the 71-7 transition form factor as well as the pion's 
elastic form factor which are calculated in the few GeV region within the modified HSA. 
This success nicely demonstrates the universality of the pion wave function, i. e., its process 
independence. Of course, the present quality of the form factor data does not pin down 
the form of the wave function precisely. Wave functions close to the asymptotic form are 
compatible with the data likewise. The BHL form is an example of such a wave function. 
On the other hand, the CZ wave function as well as other strongly end-point concentrated 
wave functions are clearly in conflict with the data and should therefore be discarded. The 
use of such wave functions in the analyses of other exclusive reactions, e. g., 77 — > nn or 
B — > 7T7T, is unjustified and likely leads to overestimates of the perturbative contributions. 
We also analyzed the 77-7 and rj'-^y transition form factors within the modified HSA. As- 
suming for the SU(3) basis states the same asymptotic wave functions as for the pion, we 
determined the decay constants and the mixing angle from a combined fit to the data on 
form factors and decay widths and found: f v = 175 MeV, f v i = 95MeV and i?p = —18°. 
These values correspond to / 8 = 145 MeV which is smaller than the chiral perturbation 
theory prediction. Keeping / 8 at the predicted value of 163 MeV, we found: f v = 203 MeV, 
f v r = 76 MeV and dp = —21°. Owing to the experimental uncertainties in the rj — > 77 
and of the moderate difference in the y 2 we cannot exclude the second set of parameters 
although the first set is favored from our analysis. 

The parameters are determined under the assumption that the asymptotic wave function is 
close to reality. Since our values of the parameters are fairly compatible with those found 
in other analyses, we are tempted to conclude that the 1]$ and rji wave functions we use, are 
indeed approximately correct. 
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FIGURES 



FIG. 1. The basic diagrams for the meson-photon transition form factor. 

FIG. 2. The n-j transition form factor vs. Q 2 . The solid (dashed) line represents the predic- 
tion obtained with the modified HSA using the asymptotic (CZ) wave function. The dotted line 
represents the results of the standard HSA (for the asymptotic wave function). Data are taken 
from [ 10 1 . 



FIG. 3. The 77-7 and rf-j transition form factors vs. Q 2 . The solid lines represent the predictions 
obtained from the modified HSA using the asymptotic wave function. Data are taken from PLUTO 
1 (■), TPC/27 § (•) and CELLO @ (o). 

FIG. 4. The elastic form factor of the pion vs. Q 2 . The dash-dotted (dash-dot-dotted) line 
represents the perturbative contributions obtained from ( |4.1| ) using the asymptotic (CZ) wave 
function. The solid (dashed) line represents the overlap contribution obtained from the soft part 
of the asymptotic (CZ) wave function. For comparison the overlap contribution obtained from the 



BHL wave function is also shown (dotted line). Data are taken from |3C 



FIG. 5. The valence quark distribution function xq^(x) at Q 2 ~ 5GeV 2 . The shadowed band 
represents the parameterization by |37j . For other symbols refer to Fig. pi The bar and the arrow 
indicate the range x in which an extrapolation for the proton structure function is used. 
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